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A dynamical model based on a continuous addition of colored shot noises is presented. The 
resulting process is colored and non-Gaussian. A general expression for the characteristic function 
of the process is obtained, which, after a scaling assumption, takes on a form that is the basis of the 
results derived in the rest of the paper. One of these is an expansion for the cumulants, which are 
all finite, subject to mild conditions on the functions defining the process. This is in contrast with 
• the Levy distribution -which can be obtained from our model in certain limits- which has no finite 

moments. The evaluation of the power spectrum and the form of the probability density function 
in the tails of the distribution shows that the model exhibits a 1// spectrum and long tails in a 
natural way. A careful analysis of the characteristic function shows that it may be separated into 
a part representing a Levy processes together with another part representing the deviation of our 
model from the Levy process. This allows our process to be viewed as a generalization of the Levy 
' process which has finite moments. 
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£ ■ I. INTRODUCTION 

o 
o 



The nature of the probability distribution of stock market prices has been discussed quantitatively for over a 
century [jj]. An early conjecture that the distribution was Gaussian was found not to be a good fit, largely because of 
the long tails found in financial data Q . A later suggestion that the Levy distribution was a better fit seemed more 
promising, since this distribution does at least have long tails ||. On the other hand this distribution has no finite 
moments, which is a severe limitation. The solution of truncating the distribution in order to obtain finite moments 
is rather ad hoc and artificial [Q. Moreover, the Levy distribution has "too fat tails" -as opposed to the "too thin 
. tails" of the Gaussian- to give a good fit to data || . 

This paper explores in more detail a model previously introduced |6) in order to resolve these difficulties and 
give a more complete explanation of the appearance of non-Gaussian and self-similar fat tails in the probability 
distribution ||, while still keeping the important feature that all moments are finite. The model has the Levy 
distribution as a limiting form and has the form of a type of "Edgeworth expansion" for the Levy distribution. We 
will give the explicit form of this relationship later in the paper (the Edgeworth series is an expansion procedure which 
gives corrections to the Gaussian distribution in those cases where the Central Limit Theorem applies Q). 

While the motivation for the model was the explanation of stock market data, we want to stress in this paper 
the more general features of the model which we expect will have applications in other areas including physics. An 
important aspect of our integrated process is that it exhibits a 1/f spectrum. We say that noise is 1/f if any 
correlated random process has a power spectral density which is, in some range, inversely proportional to a power of 
the frequency, i.e. 1/ f v (y > 0). Strictly speaking 1/f noise corresponds to the case when v ~ 1 otherwise one talks of 
"flicker noise" . However, here we use the term 1 // noise to apply to any power law with v > 0. From a historical point 
of view the first experimental observation of 1/f noise was made in 1925 by J. B. Johnson in studying non-stationary 
currents in vacuum tubes and the first attempt of a theoretical explanation was given by Schottky § in 1926. 
Since then 1/f noise has appeared, not only in semiconductor physics jlC| ], but in many areas of physics and even 
natural and social sciences [|llj- [jl4|. In spite of its ubiquity 1/f noise is not well understood from a theoretical point 
of view, and although several mechanisms to generate the noise have been proposed [jllj- |l2"| , as is pointed out in |l5| ] 
they appear to be very specialized and do not address the universality of this noise. In this sense, we hope that the 
approach herein will shed some light in the understanding oil/ f noise as it can be considered a physical mechanism 
generating such kind of noises. 

The paper is organized as follows. In Sect. II we present the main results of the paper for clarity, since many details 
are quite technical. In Sect III we describe in detail our integrated process. In Sect. IV we impose the self-scaling 
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property on the probability distribution and obtain several relevant functions. Section V is devoted to moments and 
cumulants and the evaluation of the power spectrum which shows the 1// character of our process. In Section VI we 
study the asymptotic behavior of the distribution while in Sect. VII we present, starting form our distribution, an 
Edgeworth-type series for the Levy distribution. Conclusions are drawn in Sect. VIII and some more technical details 
are in Appendices. 



II. MAIN RESULTS 



In this section we wish to state the main results of the paper without giving the derivations or being careful to 
state the range of validity for which they hold. We hope that this will give a good indication of the scope and nature 
of our results; the interested reader can then fill out this basic framework by proceeding to later sections. 

As explained in the next section, our integrated process X{t) is a continuos superposition of colored shot noises 
parametrized by u. It takes the form 



X(t) 



5^A fc (u)0(t-T fc (u);u) 



.fc=i 



du. 



There are several components: (i) The pulse shape function, (f)(t,u), (ii) the jump amplitude Ak(u) of the fcth pulse, 
(iii) the jump time Tk(u) of the kth pulse. The jump amplitudes, A/-(u), are identically distributed random variables 
distributed according the probability density function (pdf) h(x,u), i.e. h(x,u)dx = Probja; < Ak(u) < x + dx}. The 
jump times are assumed to follow a Poisson process with parameter X(u). Therefore the model is characterized by the 
functions <fi(t,u), h(x,u) and \(u). There are very few restrictions on these functions (one is the causality condition 
<j>(t, u) = for t < 0) but we will argue that it is natural to assume the scaling forms 



<f>(t,u) = <j>[\(u)t] 



and 



h(x, u) 



a(u) 



a{u) 



where <j(u) is the standard deviation of the jump amplitudes. The model is now specified by the four functions of a 
single variable <f>, h, a and A. 

One of the main goals of this paper is to find the probability distribution of the process X(t), this distribution will 
be obtained through the characteristic function (cf) p(u>,t) = (exp[iu>X(t)]), which is the Fourier transform of the 
density p(x, t). By assuming the scaling forms given above the one time distribution of X(t) is explicitly given by the 
following cf 



p(to, t) = exp i —btuj 



'"f 

Jo 







ds 




1- / h 




} 


Jo 









dz 

where b > 0. When <p(x) — 9(x) is the Heaviside step function then the input shot noises are white and X(t) is the 
Levy process: 



p(u),t) = e 



-MtLU° 



where M = b J °° dz[l — h(z)]/z 1+a . This provides us with an alternative interpretation of the Levy process since, in 
our case, Levy processes are continuos superpositions of families of white Poissonian shot noises. 

Contrary to Levy processes, our integrated process can have finite moments of any order, thus cumulants are given 

by 



«*"(*)» - 



(bt) n/a 

Jo 



<P n (x) 

x n/a 



dx, 



(n = 1,2,3,- • •). Note that the second cumulant, (i.e., the variance), is proportional to {bt) 2 / a and X(t) presents 
anomalous diffusion behavior (see Sect. V for a detailed discussion on limiting values and bounds for exponent a and 
other parameters related to the asymptotic behavior of the pulse shape function (j>{x)). 

We define the stationary correlation function by C(r) = lim 4 ^ 00 (V(t + r)X{t)). For our process this reads 



f°° dz f°° 

C ^ = b J —J <f>(bz- a t')<t>{bz- a (t' + r))dt'. 
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The power spectral density of the process X(t), given by the Fourier transform of the stationary correlation, is 



CM 



K 



and X(t) is 1/ / noise with exponent v = 1 + 2/a. 

We also can perform the asymptotic analysis of the probability distribution of X(t) without having to specify any 
particular form for h(x) and 4>(x) thus keeping the maximum level of generality. Specifically we show in Sect. VI that 
the center of the distribution is approached by a Levy distribution 



p{uj,t) 



-L(t)u 



where < S < 2. We refer the reader to see Sect. VI for more details and for the behavior of the tails of the 
distribution which are mainly determined by the behavior of the jump pdf h(x). 

The relation to the Levy distribution is explored in more detail in Sect. VII where we present and alternative (and 
exact) expression for the cf which decompose the distribution of the integrated process X(t) into that of Levy plus 
an additional term: 



lnp(w,t) = lnp Lcvy (w,i) 



1 



~h([bt/x] 1/a io) - ~h([bt/x}^ a LU(l>(x)) 



dx 



l-h^bt/xY^uj^x)) 



dx. 



Note that when 4>( x ) is the rectangular step function this equation reduces to the Levy distribution. Therefore, when 
4>(x) is a step-like function close to the Heaviside function, this alternative expression can be used as the starting 
point of an Edgeworth-type expansion procedure giving corrections to the Levy distribution. 



III. THE INTEGRATED PROCESS 



Let X(t) be a random process formed by a continuos superposition of independent shot-noise processes: 

/>oo 

X(t) = / Y(u,t)du, (1) 



where for any fixed time t, Y(u, t) are independent random variables for different values of parameter u (see Eq. ( |l(i| ) 
below) and for any fixed value of it, Y(u,t) is a colored shot-noise process represented by a countable superposition 
of pulses of identical shape: 

00 

Y(u,t)=J2Mu)cl>[t-T k (u);u], (2) 

k=l 

where Tk(u) marks the onset of the fcth pulse, and Ak(u) is its amplitude. Both Tk{u) and Ak{u) are independent and 
identically distributed random variables with probability density functions given by h(a,u) and tp(t,u), respectively. 
The pulse shape 4>(t,u) has to fulfill the "causality condition", i.e., 4>{t,u) — for t < flT6j , 

We assume that the occurrence of jumps is a Poisson process, in this case the shot-noise Y(t, u) is Markovian, and 
the pdf for the time interval between jumps, ^{t, u)dt = P{t < Tk(u) — Tk-i(u) < t + dt}, is exponential: 

V»(t,u) = \{u)e- x{u)t (i>0), (3) 

where X(u) is the mean jump frequency, i.e., 1/X(u) is the mean time between two consecutive jumps Jl7| . We recall 
that jump amplitudes Ak(u) are identically distributed (for all fc = 1, 2, 3, • • •) and independent random variables (for 
all k and u). In what follows we will assume that they have zero mean and a pdf, h(x, u)dx — Y{x < Ak (u) < x + dx}, 
depending on a single "dimensional" parameter which, without loss of generality, we assume to be the standard 
deviation of jumps cr(u) = \J [Aj,(u)). That is, 

1 [ x 1 

h(x,u) = —j—rh —— . (4) 
<j(u) l<J(u) 
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Before proceeding further with the probability distribution of the integrated process X(t) given by Eq. ([j]), we note 
that following Rice's method jl6) one can easily obtain all the probability distributions of the shot noise Y{t,u) via 
their cf's 



PY(oJi,h; ■ ■ ■ \ u n ,t n ;u) = ^exp 

In Appendix A we show that 

hxpyfa, t; u) = —X(u) 

and (supposing that ti > ii) 

\np Y (ui,ti;U2,t2;u) = -X(u) 



fe=i 



t- / h[uja(u)(j)(t',u)}dt' 



(5) 



(6) 



h[u)i<j{u)4>(t' , u) + u>2<?(u)(f)(t' + t 2 — ti,u)]dt' 

h[ui2<7(u)(j)(t\u)]dt' 



(7) 



where h(uj) is the Fourier transform of the jump pdf h(x). 

Let us now evaluate the probability distribution of the integrated process X(t). In terms of the cumulants, ((Y(t, u))) 
of the shot noise Y(t, u) we see that the one time characteristic function of X(t) can be written as 



px(u,t) = exp< ^2 



"fcT 



That is, 



,fc=i 



°° (iuj) k f°° f°° 
lnp x (u,t) = ^±-f- / .-. / ((Y(u 1 ,t)---Y(u k ,t)))du 1 ---du k 

k=l J ~oo J — OO 



Y(u, t)du 



But, by our assumptions on the process Y(u, t) we have 

«F(ui, t) ■ ■ ■ Y(u k , t))) = ((Y k { Ul ,t)))6{ Ul - u 2 ) ■ ■ ■ 5K-i - u fc ). 

Therefore, 



(8) 



(9) 



(10) 



fe=i 



fc! 



fc POO 



((Y k (u,t)))du, 



that is, 

/oo 
-oo 

Note that this line of reasoning can be easily extend to the nth time distribution, with the result: 

I' oc 

lnpx(wi,ti;- •• = / lnp Y (ui,ti;- ■ ■ ;u n ,t n ;u)du. 



(11) 



(12) 



Going back to our integrated process we have from Eqs. ©-(Q) and ([TT])-(|l2]) that the one time characteristic 
function of X{t) reads 



lnp(u),t) = — J duX(u) 
while the two time cf is (t 2 > t\) 



t- I h[u)a{u)(j){t' ,u)]dt' 



(13) 
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h) P!w'[ . ' [ : „J j. 1 2 ! = / du\{u)\ J dt' 



h[uji<j(u)<fr(t' , u) + uj2cr{u)(j){t' + ti — ti,u)] — 1 
ta-ti 



rff' 



%2<7(u)<£(t',u)] - 1 



(14) 



where we have dropped the subscript X. Obviously these are formal expressions, as long as we do not provide the 
functional dependence of A(u) and cr(tt) on the parameter u. We will do so in the next section using scaling arguments. 

IV. SCALING 

In order to proceed further we need to specify the functional form of A(u) and a(u). Of course that form will depend 
on the specific features of the problem at hand. At this point we choose what seems to us one of the most general 
ways of proceeding, we thus suppose that our integrated process X{t) possesses self-scaling properties. Following this 
path we must first assume that the pulse function is of the form 

4>(u,t) = 0[\(u)t}, (15) 

which turns <j>(u,t) into a function of the single dimensionless variable \{u)t. Substituting this into Eq. (|l~3|), defining 
new integration variables s = t' /t and z = u>a(u), and supposing that er(— oo) = and er(oo) = oo, we obtain 

]np(u,t) = - J dz ~l j 1 ~ J h[z(t>(^st)]dsX , (16) 

where the prime on a denotes derivative. We now impose the self-scaling property on the cf, that is, we assume that 
p(ui,t) is a function of the single variable Lot 1 / a : 

p(u>,t) = f^t 1 /"). (17) 

On the other hand we note that in Eq. ([if]), the quantities a' and A are functions of z and lu. Then scaling ( |l7| ) 
implies 

X = B(z)uj a , — = A(z)Lj a , (18) 
uja' 

where A(z) and B(z) are arbitrary functions to be determined. From these two relations we get a' — C(z)/u>, where 
C(z) = B(z) j 'A(z). In the original variable u we have 

a' = (19) 

But a' = cr'(u) is independent of u). Therefore, the unknown function C has to be of the form C{uj/j) = kuja where k 
is a constant. Hence <j'(u) = ka(u), whence cr(u) = <Joe ku . Finally absorbing the constant k inside the variable u we 
obtain the functional dependence of the jump variance a on parameter u: 

a(u) = cr e", (20) 

where Co is a constant. Moreover we see from Eq. ( |l8| ) that A = B{u>a{u))u> a . But, again A = A (it) is independent 
of lu. In consequence B(coa) — b (Lua)~ a where b is an arbitrary constant. Substituting this into the first relation of 
Eq. ( |l8| ) yields the "dispersion relation" between the mean frequency A and the jump variance a: 

X = b/<j a . (21) 

The functional dependence of A on u is obtained by combining Eqs. ([20|)-(pT|), 

A(u) = \ e~ au , (22) 

where Ao = &/co> or equivalently 

b = \ Q a%. (23) 
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Collecting results we see from Eq. (|lq) and Eqs. (20)-(E3|) that the one-time characteristic function reads 



p(uj, t) = exp < —btcj a 



dz 

7 l+a 



h[z(f> {btuj a s/z a )}ds 



It is sometimes convenient to rewrite this equation and use the alternative form of p(uj,t) given by 



or equivalently 



p(u), t) = exp < — bt 



p(ll>, t) = exp < b 



dz 

yl+a 



dz 



1 - / h[zuj(j){bts/z a )]ds 



dt' 



h[zuj(f>{bt'/z a )} - 1 



(24) 



(25) 



(26) 



Starting from Eq. (|l4j) and following an analogous reasoning based on the scaling assumption, we obtain the 
following expression for the two time characteristic function of the integrated process (with t% > tx) 



lnj5(wi,ii; UJ2M) = b 



f dz \ 




{!>' 


Jo ^ 1+a ] 





h[zuj 1( l>(bz- a t') + zu 2 (t>{bz- a {t' +t 2 - h))] - 1 



dt' 



h[zuj 2 (j){bz- a t')] - 1 



(27) 



Equations (|24|)-(|27|) are some of the key results of the paper. Since, as we will see next, they constitute a generalization 
of the Levy distribution with finite moments. 



V. MOMENTS, CUMULANTS AND POWER SPECTRUM 

We first note from Eq. (kJ) that if the pulse shape function is the Heaviside step function: 



1, if <>0 
0, otherwise, 



then the integrated process X{t) is identically a Levy process, regardless the jump pdf h(x): 

p(uj,t) = exp [-Mtuj a ] , 

where 

dz 



M 



1 - h{z) 



(28) 



(29) 



(30) 



Therefore, following our model, Levy processes can be viewed as a continuos superposition of families of rectangular 
pulses occurring at random Poisson times. The usual range of the exponent a in Levy flights is < a < 2. In such 
a case X(t) has no finite moment but the first one flTsfl . In actual situations, one is unlikely to meet with perfect 
rectangular pulses (showing sudden changes) in such a case all moments can be finite and are easily evaluated from 
Eq. (E9) through 



Thus for instance the second moment is given by (note that due to Eq. (|) h"(0) = -1) 



(X 2 (t)) = bt 



dz 



z 



a-l 



cj) 2 {bts/z a )ds. 



(31) 



As an illustrative example suppose that our pulse function has the form 



G 



m = \ 1 ~ e ~ kt ' ' li i >0 (32) 

w \ 0, otherwise, v ; 

where k > is a constant (note that if k is large then <fi(t) approaches to the rectangular pulse (|2§|)). In the Appendix 
B we show that (see also Eq. ( |43| ) below) 

(X 2 (t)) = Dt 2/a , (2>a>2/3), (33) 

where D = ab 2 k- 1+2 ^ a (l - 2 2 / Q ~ 2 )r(2 - 2/a)/(2 - a). Since 1 < (2/a) < 3, Eq. © clearly shows a superdiffusivc 
behavior. 

In fact we can easily obtain a closed expression not for the nth moment but for the nth cumulant 

((X™(t)))^ l -™ a " ln ^' i) 



duj n 

^From Eq. ( p5| ) we have 

/>00 />1 

((X n (t))) =r n bfh ( - n \0) / z n - 1 - a dz <t) n (Us/z a )ds. 

Jo Jo 

If in the double integral on the right hand side of this equation we define a new integration variable x by s — [z a /bt)x 
and exchange the order of integration we get 



(bt/x 



z n ~ x - a dz I 4> n (bts/z a )ds = — / cj) n {x)dx I z n ~ x dz, 
o Jo bt J J 

but the last integral is trivially evaluated, and for the nth cumulant we have 

n Jo x i 

Taking into account that ^'■"■'(O) = for n odd (we have assumed a symmetric jump distribution h(x)) we write 

((X 2n -\t)))=0, (35) 

and 

n^t^rp^^j-^ (36) 

(n = 1, 2, 3, • • •). In order to check the convergence of these expressions and the existence of moments, we first assume 
that 

<p(x) - x , (x -» 0), (37) 



(/3 > 0) then the convergence of the integral on the right hand side of Eq. (|36|) as x — > implies that the scaling 
exponent a has a lower bound: 

" > 2 " a . (38) 
1 + 2n/3 v ; 

On the other hand if we assume that 

<Xa;)~x 7 , (a: -►£»), (39) 
then the convergence of (|3^ ) when x — > oo implies that the scaling exponent a also has an upper bound: 

1 1 

-> 1+ —. (40) 

Moreover when 7 > then if 
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all cumulants will exist. Note that Eq. (|4^) holds whenever 7 > — l/2n. Therefore, 7 > is a sufficient condition for 
its validity. On the other hand if 7 < —1/2 there is no upper bound on the accepted values of a. We also observe 
that for a step-like function, as that of Eq. ( 32 ) , where 7 = then all cumulants will exist if 

2 > a > 1/(3. 

Finally, for any integrable function <j)(t) over [0, 00) there is no upper bound for a and the only condition on a for 
having all moments finite is that a > 1/(3. 

We close this discussion on moments and cumulants with an example. Suppose that the pulse function is given 
by the step-like function (|32|). In this case 7 = 0, (3 = 1 and all moments (and cumulants) will exist if 1 < a < 2. 
Cumulants are given by Eqs. (B5|)-(M). In Appendix B we show that 



( 1 _ p -kx\n 

■i '—da = A„fc- 1+ ™/ Q r(l + n - n/a), 



(42) 



where the numbers A n are given by Eq. (B3) of Appendix B. Finally, 



{{X 2n (t))) = i '— — y -^A 2n k- 1+2n/o T(l + 2n- 2n/a){btf n l a . (43) 

We finish this section evaluating the power spectrum of the integrated process X(t). Let us first evaluate the 
correlation function 



d 2 

(X(t + r)X(t)) = - p( Wl) i; U2, t + t) 

OUJ1UJ2 



U)\ =io'2 =0 



,/From Eq. ( P7| ) we get 

(X(t + r)X(t)) =b -JL <f>(bz- a t')cl>(bz- a (t> + r)). (44) 
Jo z Jo 

Let C(t) be the correlation function in the stationary limit t — > 00, i.e. 



^From Eq. (Q) we have 



C(r) = lim (X{t + r)X(t)). 

t — fOO 



= b I i / <Kbz~ a lf)tKbz- a y + r)). (45) 



Note that the (stationary) variance C(0) = 00 which agrees with the superdiffusivc behavior of X(t) given by Eq. (|3^). 

The power spectral density of our process is thus given by the Fourier transform of the stationary correlation 
function 

00 

ILJT / 



C{U)) = / e-' lUT C(T)dT. 

J —00 

Substituting Eq. (|4|) into this equation, performing simple changes of variables and taking into account the causality 
of the pulse function <p(t) we finally obtain 

<?M = (46) 

where 

,2/a roo 

K=— e /a \m\ 2 d^ (47) 
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and 0(£) is the Fourier transform of (pit). The power spectral density C(uS) exists if the integral 

j= [°° e /a \m 2 dt < oo. 

Jo 

In order to prove the existence of J we first need that the Fourier transform of the pulse function, exists. Note 

that any step or step-like function does not have a Fourier transform and consequently the power spectrum is infinite. 
For the existence of (/>(£) it suffices that <p(x) be absolutely integrable, and from the asymptotic behavior given by 
Eq. @: 

4>{x) ~ x 1 [x — * oo), 

we have to impose that 7 < — 1. This in turn implies that </>(£) ~ £~ 1_ 7 as £ — > and, since 7 < — 1, the integral J at 
its lower limit is always finite for any a > 0. On the other hand if (f>(x) satisfies Eq. (p7\) as x — * 0, then 0(£) ~ t^ 1 -! 3 



as £ — * 00. Hence, J will be finite if a > l/(/3+ 1/2) (see also Eq. (38)). Therefore, the process X(t) has a finite 



power spectrum if (fs(x) is absolutely integrable on the real line, and X(t) has a finite second cumulant, Eq. (38), i.e 



a > /JTT72- ^ 

In such a case, we see from Eq. (|46| ) that -X~(t) is 1// noise with exponent v = 1 + 2/a. Moreover (recall that a > 0, 
and Eq. @) 

1 < v < 2(1 + j8). (49) 

Finally, when <fi(x) is analytic at x = then /? = 1, 2, 3, • • • is a positive integer and X(t) is flicker noise l/f v with 
1 < z/< 2(1 + n) (n = 1,2,3,- • •). 

VI. ASYMPTOTIC BEHAVIOR 

We will now examine the asymptotic behavior of the one time probability density function (pdf) of the integrated 
process X(t), p(x,t)dx = P{x < X(t) < x + dx). For this analysis we distinguish two regions: the "center" (x — > 0) 
and the "tails" (x — » ±00) of the distribution. We cannot have a closed expression for the pdf p(x, t) until the pulse 
function <p(x) and the jump pdf h{x) are both specified. Therefore, we will perform the asymptotics on the cf p(lo, t). 
As a well known feature of the harmonic analysis the center of the distribution is determined by the large us behavior 
of the cf, while the tails are determined by p(u>, t) when us — > p9[ |. 

We deal first with the center of the distribution where us — > 00. If we assume that the pulse function 4>(x), as 
x — > 00, satisfies Eq. (^9h: 



then 



4>(x) ~ X 7 , [x — > 00), 



h [z4> (btus a s/z a )] - h [(btus a y z 1 - a ' f s' y ] , (us 



Substituting this into Eq. (G4) and performing the change of variables £ = (btus") 1 z 1 ai we obtain the following Levy 
distribution: 

p(us,t)^exp{~L(bt) 1 ^ 1 - a ^us a /^ a ^Y (us^oc), (50) 



where 



L = -J— / — [l - / h{^)ds 

l-a7 JO £( 1 +«-"7)/(l-"7) [ Jo 



(51) 



Note that due to the bounds discussed above (see Eq. (|40|)) we have 1 — 07 > a/2 hence the Levy exponent in Eq. (|5_ 
satisfies 

< — 5L_ < 2 , (52) 
1 — 07 
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and Eq. (|5(j) is well denned. We also note that for a step-like pulse function <fi(t) where 7 = we obtain the same 
Levy distribution, Eq. ((2^), that satisfies the model for sudden pulses (^8|). Therefore, for any pulse shape function 
satisfying condition ( |39| ) the center of the pdf is given by a Levy distribution. 

Let us now obtain an asymptotic expression of the pdf p{x,t) when x — > ±00, which will be valid if 4>(t) obeys 
Eqs. (p7]) and (p39), and the exponent a is bounded by 



/T 1 < a < 1/(7 + 1/2). 

In this case, we see from Eq. ( fll|) that all cumulants exist. So taking the u> — ► limit of Eq (]25| ) we find 

dz 



p(u), t) ~ 1 - bt 
The Fourier inversion of Eq. (|5^) yields 

p(a;, <) ~ bt 



1 - 



h[zu)(/){bts/z a )}dz 



dz 

yl+a 



h(x,z,s)dz (x — > ±00), 



(53) 



(54) 



where we have dropped delta function terms which have no contribution as x — > 00. Moreover assuming symmetric 
jump distributions h{x) we have 



1 



h(x,z,s) = — / h [zuxfi (bts/z a )] cos ux dw 



1 



z(j){bts/z a ) 



z(j)(bts/z a ) 



Substituting this into Eq. (b3) we see that 



p(x, t) ~ bt 



dz 

y 2+a 



ds 



o z 2 +«y (t>{bts/z<*) 



z(/)(bts/z a ) 



(x —>■ ±00). 



(55) 



Therefore the tails of the distribution are determined by the jump pdf h(x) and the pulse shape function 4>{t). 
Finally for the rectangular pulse (Eq) we have 



p(x,t) 



bt 



ll + Q 



y a h(y)dy, 



(56) 



which agrees with the expected tail behavior of the Levy distribution pi 



VII. RELATION TO THE LEVY DISTRIBUTION 

In section V we obtained general expressions (|3^)-(|3^) for all of the cumulants, from which it follows that 

lnp( W ,t) = X; ! - r «^ n (*)» 

71— \ 

^ 2n(2ri)\ J x 2n / a K ' 

In addition, in section VI, wc showed that p(x, t) was a Levy distribution at the center of the distribution (x — » 0) 
and took the form (|5^) in the tails (x — > ±00) of the distribution. In this section we will show how the distribution 
can be separated into a Levy distribution plus an additional term. This term takes the form of a single integral which 
can be evaluated once the functions <j) an d h have been specified. 

We begin the analysis by changing variables from s to x = bts/z a (z fixed) in Eq. (p5[). This gives 

r°° dz f bt/z 



f°° dz f 0l/z r - 1 
hip(u), t) = — I — / dx \ l — h \zu>(j){x)} \ 
Jo z Jo L J 



^ 1/ 

dx 
[i Jo 



(b * /x)1/ ° dz ( ~ i 

— {l-h[zu)4>(x)]j (58) 
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changing the order of integration. At this point we factor out the contribution from the Levy process by writing (|56 
as 



dx 



dz r - i f°° r( bt / x ) 1/a dz f- ~ i 

— |l — h [zlu] | — J dx J — | h [zlu] — h [zu<f>(x)] j 



or, after defining 



g{w,t) = 1 - h(u,t), 



(59) 



as 



dx 



(bt/x) 1/o 



dz 



oo j-ibt/x) 1 '" dz 

dx 



{<J [zu] -g[zu<f>(x)]}. 



The first term is just lnp(w, t) for the Levy processes (see Eqs. (P9j)-(|30|))- The second term can be simplified by first 
writing it as 



dx 



(bt/x 



){zuj) 



(bt/x) 1 /"^ z 



dz 



and then integrating by parts to give 



g(zui) 



dz 



dxx^— 



g(zuj) 



dx J(bt/xy/^4> z 



dz. 



I(bt/x) 1 /°> l p z 

We assume that h(u>) is analytic at li — and integrable, then 

g(ui) ~ id 2 (ui — > 0) and <?( w ) ~~ * 1 (w — ► oo), 
and since < a < 2 the first term in Eq. (BG) is zero. Finally 



(60) 



1 



lnp(w, t) - lnp Lcvy ( W , t) + - / g([bt/x]^ a oj) - g([bt/x]V a u<l>(x)) 



dx 



+ [™ x^rlHlbt/x} 1/a ^(x))dx, 
Jo 9\ x ) 



(61) 



where cj)'(x) is the derivative of the pulse shape function and 

lnpLevy(w,i) = -Mtw" 

where M is given by Eq. (^). Note that when (j>{x) is the Heaviside step function the integrals on the right hand 
side of Eq. ( pH ) vanish and Eq. (|6l] ) reduces to the Levy distribution. Therefore, we can look at the second term on 
the right hand side of Eq. (^TJ) as a correction to the Levy distribution when <fi(x) is not exactly a Heaviside function 
but a step-like function very close to the Heaviside function. This may be evaluated, in principle, for any given <p and 
h. For instance we could take cf> to be of the form (|3^) with k large and the Lorentzian 

u; 2 /2 



l + c^ 2 /2' 



corresponding to h(x) — e *^\ x \ /y/2. The form of the correction terms depends on the choice of the functions to an 
extent, and so we will not discuss the explicit form it takes here. 



VIII. CONCLUSIONS 



In this paper we have presented and analyzed a dynamical model based on a process which is a superposition of 
colored Poisson noises. The model was shown to have several attractive features. The probability density function has 
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long tails which emerged in a natural way and, unlike the Levy distribution, all the moments of the distribution are 
finite. We believe that these properties make the distribution an ideal candidate for describing stock market prices [0 . 

A property what may have more relevance to physics and other natural sciences is the appearance of 1// noise in 
the power spectrum. Once again we would like to stress that this result flowed naturally from the nature of the model 
and the scaling assumptions which reduce h and 4> from functions of two variables to functions af a single variable. 

In a more mathematical context, we believe that the decomposition of the cf of our model into that of the Levy 
plus additional terms is interesting, both as an example of an Edgeworth-type expansion and for the nature of the 
corrections to the Levy distributions when the parameters of our model are chosen so that our distribution is near to 
the Levy one. 

There are still some open questions. One of them is the extension of the model to the increments of the process 
Z(t, t — to) = X(t — t + t) — X(t — to) (t > to), since in this case we believe that the process Z(t, t — to) becomes 
stationary when it starts in the infinite past (to — > oo). Another interesting and open question is the actual application 
of the model to financial time series where some non- white correlation is observed p0| |. Both points are presently 
being investigated. 
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APPENDIX A: CHARACTERISTIC FUNCTION FOR COLORED SHOT NOISE 

By generalizing Rice's method we will now obtain the probability distribution of the shot noise Y{u, t) defined 
by Eq. (§): 

OO 

Y(u,t)=J2 A k(u)<f>[t-T k (u);u], (Al) 
fc=i 

where we assume that the random variables A^(u) and T^(u) are identically distributed and statistically independent. 
The jump amplitudes are described by the pdf h(x,u)dx = Probja; < Ak(u) < x + dx} and the jump times Tk(u) 
follow a Poisson distribution of parameter X(u). Define 

p(xi,ti;x2, t 2 ; u)dx\dx 2 = Probjzi < Y(u, ti) < x\ + dx\;x 2 < Y(u, t 2 ) < x 2 + dx 2 } 

to be the joint pdf of the process with t 2 >t\. This pdf can be written as 



oo oo 



p(x 1 ,t 1 ;x 2 ,t 2 ;u) = ^ ^2 P(xi,t 1 ;x 2 ,t 2 ;u\n 1 ,n 2 )P(n 1 ,t 1 ;n 2 ,t 2 ;u), (A2) 

ni=0 ri2— 

where p(x\, t\\ x 2 , t 2 ; u\n\, n 2 ) is the conditional pdf assuming that exactly n\ pulses have occurred at time t\ and 
n 2 pulses at time t 2 . P(ni,ti;n 2 ,t 2 ;u) is the joint probability for the occurrence of such pulses. Since t 2 > t\ then 
n 2 > n\ and 

P(ni,t 1 ;n 2 ,t 2 ;«) = /^ n2 - ni; * 2 - <i;u ) P ( ni '* i;tt )' if " x ~. n2 (A3) 
10, otherwise, 



where 



P(m,T-,u)= M u W m e -Hu)r (A4) 



is the Poisson distribution. Substituting Eq. ( |A3| ) into Eq. (A2) and defining t\ — t, t 2 = t + At, n\ = n and 
n 2 — m = to, we obtain 

oo oo 

p(x\, t;x 2 ,t + At; u) = p(x\, t; x 2 , t + At; u\n, n + m)P(m, At; u)P(n, t; u), 

n—0 m— 
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and the characteristic function reads 

oo oo 

p(u>\, t; u>2, t + At; u) = p(w%, t; u>2,t + At; u\n, n + m)P{m, At; u)P(n, t; u). (A5) 

n— m— 

Note that p(u>i, t; u>2, t + At; u\n, n + m) is the joint characteristic function of the truncated process 



Y n {u,t) =^A k (u)4>[t-T k {u);u] 



fc=i 



Hence 



p(u>i,t; ui 2 ,t + At; u\n, n + m) = ( exp ioj 2 ^ A k {u)4>(t + At — T k (u);u) + iu\ ^ Ai(u)(f)(t — Ti(i 



fe=i 



1=1 



n+m 



+ iuj 2 A k {u)(j)(t + At-T k {u);i 

k =n+l 

Taking into account that A k {u) and T k {u) are independent and identically distributed random variables, we have 

-i n 

p(b3i_,t',tjJ2,t + At; u\n,n + m) = (exp{iui 2 A k (u)(t)(t + At - T k (u);u) + iu>iA k (u)(/)(t - T k (u);u)}) 

-i m 

x \(exp{iuj 2 A k (u)cj)(t + At~ T k (u); u)})J , 

and, since the random times T k (u) are Poissonian, 

(exp{iu; 2 A k (u)(j)(t + At - T k (u); u) + iu) X A k {u)(j){t - T k (u);u)}) 

f* dt' 

h(a,u)da / — exp{ia[u 2 (f>(t + At — t; u) + LOi<j){t — t'; u)]} 

.In t 



1 /"'- 

- / h[u 2 (j)(t 
t Jo 



+ At - t'; u) + uji4)(t - t'; u); u] dt' , 



where h{uo; u) is the Fourier transform of the jump pdf h(x; u). Analogously 



1 



(e^p{iu 2 A k (u)(f)(t + At~T k (u);u)}) = — J h [aj 2 (/)(t + At - t! ; u); u] dt' . 



t+At 



Therefore, 



where 



p(u>i, t; cj 2 , t + At; u\n, n + m) 



F(u)2,t + At;wi,f; u) 



At 



G(uj 2 ,t + At;t;u) 



F(u 2 ,t + At;uj 1 ,t;u) = / h[uj 2 (j){t + At - t'; u) + wi<j>(t - t'; u); u] dt' , 



and 



t+At 



G(cj 2 ,t + At;t;u) = I h[oj 2 (j){t + At - t' ';u);u]dt' '. 



(A6) 



(A7) 



(A8) 



Substituting Eq. (|A6|) into Eq. (|A5| ) yields 



p(ui,t; uj 2 ,t + At;u) = 



~F(wi,t;cj2,t + At;u) 



P(n,t;u)\\ 



—G{uj 2 ,t + At;t;u) 



P(m, At; \ 
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Introducing Eq. (A4) into this and performing the resulting sums we get 

p(ui 1 ,t;u! 2 ,t + At;u) = exp{-A(w) [t - F(u 2 ,t + At;u>ut;u)] - X(u) [At - G(w 2 , t + At; w x , t; u)]} . 

Finally 



p(uj l7 t;u 2 ,t + At;u) = exp< X(u) / di' h[uJi^(t', u) + uj 2 (f>(t' + At, u)] - 1 



Ai 



+ X(u) / di' ft[w 2 0(i',u)] - 1 



(A9) 



which agrees with Eq. (0). If in @ we set cji = cj and lo 2 = 0, we obtain the one time characteristic function (||): 



p(uj, t; u) = exp < A(u) / di' ft[w</>(i', «)] — 1 



(A10) 



APPENDIX B: CUMULANTS FOR A STEP-LIKE FUNCTION 



We will derive closed expressions for the cumulants (36) 



((X 2n (t))) = ( l) n ^ (2n) (0) (bf) 2n/a 



2/7 



°° 2 "(a;) 

2n/a 



dx. 



o ^ 



when 4>(x) is the step-like function (133) : 



l-e~ kt , if t > 



0. 



otherwise. 



The substitution of Eq. (B2) into Eq. (Bl) leads us to evaluate the following integral 



(-1 _ -kx\2n 
hn = I ^ f , ' dx 



x 

00 dx 



2n/a 



x 



k 2n I dui -I du 2n I x 



2n/ a 



kx I e~ kxu du 



1 />oo 



2d 



JO 



2«(l-l/a) 



exp< 



-kx I ] >dx. 



(Bl) 



(B2) 



Define the new integration variable 



then 



-1 .1 / 2n 

J 2 „ = fc- 1+2 "/« jf d Ul • • • J du 2n 



; ,„ v -l-2n(l-l/a) 



But 



/ e -^ 2n(1 - 1/a) de = r(l + 2n-2?i/a). 
Jo 



Defining the numbers 



14 



-1 ,.\ I In 

Mn = J da X - J du 2n ( ^2 



:■„ \ -l-2n(l-l/a) 

Ui 
'"-I / 



we finally have 



C— 1 ^ n /)( 2 ™) fni 

((X 2n (i))) = i '- y ^lA 2n k- 1+2n/a T(l + 2n - 2n/a)(bt) 2n/a , 



2n 



(B3) 



(B4) 



which is Eq. (43) 
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